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ABSTRACT 

We study the excitation of density and bending waves and the associated angular mo- 
mentum transfer in gaseous disks with finite thickness by a rotating external potential. 
The disk is assumed to be isothermal in the vertical direction and has no self-gravity. 
The disk perturbations are decomposed into different modes, each characterized by the 
azimuthal index m and the vertical index n, which specifies the nodal number of the 
density perturbation along the disk normal direction. The n = modes correspond to 
the two-dimensional density waves previously studied by Goldreich & Tremaine and 
others. In a three-dimensional disk, waves can be excited at both Lindblad resonances 
(for modes with n = 0, 1, 2, • ■ •) and vertical resonances (for the n > 1 modes only). 
The torque on the disk is positive for waves excited at outer Linblad/vertical reso- 
nances and negative at inner Lindblad/vertical resonances. While the n = modes 
are evanescent around corotation, the n > 1 modes can propagate into the corotation 
region where they are damped and deposit their angular momenta. We have derived 
analytical expressions for the amplitudes of different wave modes excited at Lindblad 
and/or vertical resonances and the resulting torques on the disk. It is found that for 
n > 1, angular momentum transfer through vertical resonances is much more efficient 
than Lindblad resonances. This implies that in some situations (e.g., a circumstel- 
lar disk perturbed by a planet in an inclined orbit), vertical resonances may be an 
important channel of angular momentum transfer between the disk and the external 
potential. We have also derived new formulae for the angular momentum deposition 
at corotation and studied wave excitations at disk boundaries. 

Key words: accretion, accretion discs - hydrodynamics - waves - binaries: general - 
planetary systems 



1 INTRODUCTION 

Gravitational interaction between a gaseous disk and an an external body plays an important role in many astrophysical 
systems, including protoplanetary disks, binary stars and spiral galaxies. The external potential generates disturbances in the 
disk, shaping the structure and evolution of the disk, and these, in turn, influence the dynamics of the external object itself. 
In a classic paper, Goldreich & Tremaine (1979; hereafter GT) studied density wave excitation by an external potential in 
a two-dimensional disk and gave the formulation of angular momentum transport rate or torque at Lindblad and corotation 
resonances for disks with or without self-gravity (see also Goldreich & Tremaine 1978; Lin & Papaloizou 1979). Since then, 
numerous extensions and applications of their theory have appeared in the literature. For example, Shu, Yuan & Lissauer 
(1985) studied the effect of nonlinearities of the waves on the resonant torque. Meyer- Vernet & Sicardy (1987) examined the 
transfer of angular momentum in a disk subjected to perturbations at Lindblad resonance under various physical conditions. 
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Artymowicz (1993) derived a generalized torque formula which is useful for large azimuthal numbers. The saturation of the 
corotation resonance were investigated in detail by many authors (e.g., Balmforth & Korycansky 2001; Ogilvie & Lubow 
2003). Korycansky & Pollack (1993) performed numerical calculations of the torques. Applications of the GT theory were 
mostly focused on disk-satellite interactions, including the eccentricity and inclination evolution of the satellite's orbit (e.g., 
Goldreich & Tremaine 1980; Ward 1988; Goldreich & Sari 2003) and protoplanet migration in the Solar Nebula (e.g., Ward 
1986, 1997; Masset & Papaloizou 2003). 

A number of studies have been devoted to the three-dimensional (3D) responses of a disk to an external potential. Lubow 
(1981) analyzed wave generation by tidal force at the vertical resonance in an isothermal accretion disk and investigated 
the possibility that the resonantly driven waves can maintain a self-sustained accretion. Ogilvie (2002) generalized Lubow's 
analysis to non-isothermal disks and also considered nonlinear effects. Lubow & Pringle (1993) studied the propagation 
property of 3D axisymmetric waves in disks, and Bate et al. (2002) examined the excitation, propagation and dissipation of 
axisymmetric waves. Wave excitation at Lindblad resonances in thermally stratified disks was investigated by Lubow & Ogilvie 
(1998) using shearing sheet model. Takeuchi & Miyama (1998) studied wave generation at vertical resonances in isothermal 
disks by external gravity. Tanaka, Takeuchi & Ward (2002) investigated the corotation and Lindblad torque brought forth 
from the 3D interaction between a planet and an isothermal gaseous disk. The excitation of bending waves was studied by 
Papaloizou & Lin (1995) and Terquem (1998) in disks without resonances. Artymowicz (1994), Ward & Hahn (2003) and 
Tanaka & Ward (2004) investigated many aspects of resonance-driven bending waves. 

Despite all these fruitful researches, a unified description and analysis of wave excitation in 3D disks by external potentials 
are still desirable. This is the purpose of our paper. As a first step, we only consider linear theory. Our treatment allows for 
Lindblad, corotation and vertical resonances to be described within the same theoretical framework, and in the meantime, 
density waves and bending waves to be handled in an unified manner. By taking advantage of Fourier-Hermite expansion, 
different modes of perturbation for locally isothermal disks are well organized and a second-order differential equation for 
individual mode is attained (see Tanaka, Takeuchi & Ward 2002). In order to treat it mathematically, the derived equation 
which appears monstrous is pruned under different situations for those modes with the highest order of magnitude. Then, 
following the standard technique used by Goldreich & Tremaine (1979), the simplified equations are solved and analytic 
expressions for the waves excited at various locations and associated angular momentum transfer rates are calculated. 

Our paper is organized as follows. The derivation of the basic equations describing the response of a disk to an external 
potential are briefly presented in §2. The assumptions underlying our formulation are also specified there. In §3 we examine 
the dispersion relation for free waves and various resonances which may exist in the disk. Wave modes are organized according 
to the azimuthal index m and the vertical index n, with n — corresponding to the modes in a 2D disk. In §4 we study 
Lindblad resonances. §4.1 deals with the n / 1 cases, where the solutions pertaining to the waves excited at the resonances 
are found and the angular momentum transports by such waves are calculated. The n = 1 case is treated separately in §4.2, 
because in this case the Lindblad resonance coincides with the vertical resonance for a Keplerian disk. In §5 we study wave 
excitation and angular momentum transport at vertical resonances for general n. Corotation resonances are examined in §6, 
where it is necessary to treat the n > 1 and n = cases separately. In §7 we study wave excitation at disk boundaries. We 
discuss the effects of various assumptions adopted in our treatment in §8 and summarize our main result in §9. Appendix A 
contains general solutions of disk perturbations away from resonances, and Appendix B gives a specific example of angular 
momentum flow in the disk: wave excitation at a Lindblad or vertical resonance, followed by wave propagation, and finally 
wave damping at corotation. 



2 BASIC EQUATIONS 

We consider a geometrically thin gas disk and adopt cylindrical coordinates (r,0,z). The unperturbed disk has velocity 
vo = (0, rfi, 0), where the angular velocity fi = O(r) is taken to be a function of r alone. The disk is assumed to be isothermal 
in the vertical direction and non-self-gravitating. The equation of hydrostatic equilibrium (for z<r) reads 

Thus the vertical density profile is given by 

p {r,z)= ° cxp(-£ 2 /2), with Z = z/h (2) 
\/2irh 



where h = h(r) — c/Q± is the disk scale height, c = c(r) = \fpoJpo is the isothermal sound speed, a — a(r) — j dz po is the 
surface density, and fl± is the vertical oscillation frequency of the disk. 

We now consider perturbation of the disk driven by an external potential <j>. The linear perturbation equations read 

^ + (v ■ V)u + (u ■ V)v = -— VSP + % 
dt po pi 



— + (v ■ V)u + (u ■ V)v = -— X7SP + -£Vp - V0, (3) 
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^ + V ■ (p u + v 5p) 



0, 



(4) 



where Sp, 8P and u = 5v are the (Eulerian) perturbations of density, pressure and velocity, respectively. Without loss of 
generality, each perturbation variable X and the external potential <j> are assumed to have the form of a normal mode in 6 
and t 



X(r, 8, z, t) — X(r, z) exp(im# — icot), 



(5) 



where m is a nonnegative integer and consequently, lo is allowed to be either positive or negative, corresponding to the prograde 
or retrograde wave, respectively. Note that only the real part of the perturbation has physical meaning and we will not write 
out explicitly the dependence on m and ui for the amplitude X(r, z) and other related quantities in this paper. The equation 
for adiabatic perturbations is dP/dt = c 2 a dp/dt, where c s is the adiabatic sound speed. This yields 



—iui{8P — c a 8p) 
where 

u! = u> — mfi 



c a p u ■ A, 



(6) 



(7) 



(8) 



is the "Doppler-shifted" frequency, and 

_ Vpo _ Vpo 

Po clpo 

is the Schwarzschild discriminant vector. In general, \A r \ ~ 1/r may be neglected compared to \A Z \ ~ 1/h for thin disks 
(h r). In the following, we will also assume A z = 0, i.e., the disk is neutrally stratified in the vertical direction. This amounts 
to assuming c a = c (i.e., the perturbations are assumed isothermal). Equation (6) then becomes 5P = c 2 Sp. Introducing the 
enthalpy perturbation 

V = Sp/po, (9) 
Equations (3) and (4) reduce to ' 
d 



-ILUUr 



2Qug = 
2 



dr 



-ILOU0 + —U r = (T) + (f>), 

Z\l r 



• - Po 19, 
-iu—ri + - — (rp Ur 
cr r or 



im d , . 

H p u e + tt- (poUz) 

r oz 



Here k is the epicyclic frequency, defined by 
2 20. d 



r dr 



(r'fl). 



(10) 

(11) 

(12) 
(13) 

(14) 



In this paper we will consider cold, (Newtonian) Keplerian disks, for which the three characteristic frequencies, O, Cl± and 
k, are identical and equal to the Keplerian frequency Q,k = (GM/r 3 ) 1 ^ 2 . However, we continue to use different notations 
(SI, Q±, n) for them in our treatment below when possible so that the physical origins of various terms are clear. 

Following Tanaka et al. (2002) and Takeuchi & Miyama (1998) (see also Okazaki et al. 1987; Kato 2001), we expand the 
perturbations with Hermite polynomials H n : 



<j>(r,z) 




<t>n{r) 




r/(r,z) 
u r (r,z) 


- E 


Vn(r) 
u rn (r) 


Hn(Z), 


ue(r,z) 


n 


U0 n {r) 





u z (r,z) = 



u zn (r)H' n {Z), 



(15) 



where H' n = dH n /dZ. The Hermite polynomials H n (Z) 
K - ZH' n + nH n = 0, 
H'n = nH n -i, 



(-l)"e z /2 d"(e 



-Z 2 /2 



)/dZ n satisfy the following equations: 



(16) 
(17) 



t On the right-hand-side of eq. (10), wc have dropped the term 2{d\nc/dr)r). In effect, we assume c is constant in radius. Relaxing this 
assumption docs not change the results of the paper. See §8 for a discussion. 
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ZH n = H n +1 + nH n -i, 

exp{-Z 2 /2)H n H„,dZ = v^n! 5 n 



(18) 
(19) 

We note that any complete set of functions can be used as the basis of the expansion. For an example, Lubow(1981) used 
Taylor series in the expansion of vertical variable. However, choosing the Hermite polynomials as the basis set, as we shall see, 
greatly simplifies the mathematics involved, since they are eigenmodes in variable z for locally isothermal disks with a constant 
scale height h and quasi-eigenmodes for disks with a small radial variation of h. Note that since Hi = Z, Hi — Z 2 — 1, the 
n = 1 mode coincides with the bending mode studied by Papaloizou & Lin (1995) (who considered disks with no resonance), 
and the n — 2 mode is similar to the mode studied by Lubow (1981). 
With the expansion in (15), the fluid equations (13) become 



„_ d nu (n + l)(n + 2)u 

— lUJU rn — 2\lUe n = — T l»n H W n H W n +2, 

dr r r 



-iibuen + 



-ILUUzn 



c 

where 
w n =Vn + ' 
and 



-w n /h, 

d , nu, 

— In ra H 

dr r 



Urn ~t~ 



d 

2 + -j-Urn + 

dr 



-Ug r , 



o, 



(1 ■: 



d\nh 
dlnr ' 



Eliminating ug n and u zn from eqs. (20)-(23), we have 

dw n 2mQ. 
dr rCb 

dUrn 



^-iu rn + - \nw n + (n+l)(n + 2)w„ +2 ], 
lo r 



dr 



dln(rcr) mc 



dr 



2rVtCj 



rn 



n 



Urn + — [ — + I Wn + _ —(riUrn + U r ,n-2), 



(20) 

(21) 
(22) 
(23) 

(24) 

(25) 

(26) 
(27) 



where we have defined 



D ee k 2 - Co 2 = k - {to - mQ.) 2 . 



(28) 



Finally, we eliminate u rn to obtain a single second-order differential equation for rj n [see eq. (21) of Tanaka et al. 2002] t; 

2mQ / d n«r\ m 2 D(d) 2 - nSl\ 
ruj \dr D J r 2 c 2 w 2 

D AmQ,\ 



dr 2 



d . ra\ d 
— In — — J — 

dr D J dr 



r 



(— 

\dr 



2mQ.\ 
— — w n 

ruj J 



-2 + 



n(— ln- 

\dr a 



-) 



- n + l)(n + 2 — - — In — + — — ) w n+2 
\dr dr a no J 



(n — 2)w„_2 + n(2n — l)w n + n(n + l)(n + 2)w n ^ 



(29) 



Obviously, for /i = 0, different n-modes are decoupled. But even for fj, =fc 0, approximate mode separation can still be 
achieved: When we focus on a particular n-mode, its excitation at the resonance is decoupled from the other n-modes (see 
Sects. 4-6), provided that the orders of magnitudes of rjn±2 and their derivatives are not much larger than r\ n and dr/n/dr — 
we shall adopt this assumption in remainder of this paper. Note that if rj„±2 is much larger than rj n , the coupling terms must 
be kept. In this case, the problem can be treated in a sequential way. After arranging the potential cf)„'s in the order of their 
magnitudes, from large to small, we first treat the mode with the largest n ; in solving this "dominant" mode, the coupling 
terms involving the other "secondary" n-modes can be neglected. We then proceed to solve the next mode in the sequence: 
Here the coupling terms with the dominant mode may not be neglected, but since the dominant mode is already known, these 



t On the left-hand-side of eq. (29), we have dropped the term — r 1 (dfi/dr) [nw n + (n + l)(n + 2)w n +2\. This term docs not change the 
results of this paper. Also, for a Keplerian disk with c =constant, h oc r 3 / 2 and dfi/dr = 0. See §8 for a discussion. 
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coupling terms simply act as a "source" for the secondary mode and the torque formulae derived in the following sections 
(§§4-7) can be easily modified to account for the additional source. 

In the absence of self-gravity, waves excited by external potential carry angular momentum by advective transport. The 
time averaged transfer rate of the z-component of angular momentum across a cylinder of radius r is given by (see Lynden-Bell 
& Kalnajs 1972; GT; Tanaka et al. 2002) 



271 



F{r) = (r 2 j dz d8po(r,z)u r (r,9,z,t)ue(r,6,z,t)J. (30) 

Note that a positive (negative) F(r) implies angular momentum transfer outward (inward) in radius. Using u r (r,6, z,t) = 
Re[u rn H n (Z)e i(me - LJt) ], u g (r,9,z,t) = Re [ue n H n (Z)e i{me ~ ut) ], and eqs. (2) and (19), we find that the angular momentum 
flux associated with the (n, m)-mode is 

F n (r) = n\nr 2 aRe(u rn u* 0n ) (31) 

(recall that we do not explicitly write out the dependence on m). Using eqs. (20) and (21), this reduces to (see Tanaka et 
al. 2002) 



Mr) = f) 



Wn— : (n + l)(n + 2 -w„w„ +2 

dr r 



(32) 



where w n = rj n + 4> n - To simplify eq. (32) further, we shall carry out local averaging of F n (r) over a a scale much larger than 
the local wavelength |fc| _1 (see GT). As we see in the next sections, the perturbation r\ n generated by the external potential 
0o generally consists of a nonwave part rj n and a wave part fj n ; thus w n = 4> n + fj n + fj n . The cross term between (cj> n + fj n ) 
and fj n in eq. (32) can be neglected after the averaging. The coupling term (oc w n w* +2 ) between different modes can also be 
dropped because of the radial averaging and |ui n+2 | & \w n \ (see above). Thus only the wave- like perturbation carries angular 
momentum, and eq. (32) simplifies to 

K(r)w n!^ Im / dgy (33) 



D y" dr / 

In §§4-6, we will use eq. (32) or (33) to calculate the angular momentum transfer by waves excited or dissipated at various 
resonances. 



3 DISPERSION RELATION AND RESONANCES 

Before proceeding to study wave excitations, it is useful to consider local free wave solution of the form 



rj n oc exp 



/' 



k(s)ds 



(34) 



For \kr\ 3> 1, from the eq. (29), in the absence of the external potential, we find (see Okazaki et al. 1987; Kato 2001) 

{Co 2 — n 2 ){Cj 2 — nQ. 2 ± )/uj 2 = k 2 c 2 , (35) 

where we have used h — c/Q± <g r (thin disk), and m,n<g.r/h — we will be concerned with such m, n throughout this paper. 
Obviously, for n = we recover the standard density-wave dispersion relation for 2D disks without self-gravity. In Appendix 
A, we discuss general WKB solutions of eq. (29). 

At this point it is useful to consider the special resonant locations in the disk. These can be recognized by investigating 
the singular points and turning points of eq. (29) or by examining the characteristics of the dispersion relation (35). For u> > 0, 
the resonant radii are 

(i) Lindblad resonances (LRs), where D = or Co 2 — n 2 , including outer Lindblad resonance (OLR) at Co = k and inner 
Lindblad resonance (ILR) at u> = —n. The LRs arc apparent singularities of eq. (29) — we can see this from (26) and (27) 
that the wave equations are well-behaved and all physical quantities are finite at D — 0. The LRs are turning points at which 
wave trains are reflected or generated. Note that the ILR exists only for m > 2. 

(ii) Corotation resonance (CR), where Co = 0. In general, the CR is a (regular) singular point of eq. (29), except in the 
special case of n = and d(Qa/K 2 )/dr — at corotation. Some physical quantities (e.g., azimuthal velocity perturbation) 
are divergent at corotation. Physically, this sigularity signifies that a steady emission or absorption of wave action may occur 
there. Note that no CR exists for m = 0. 

(iii) Vertical resonances (VRs), where Co 2 — n£l\ (with n > 1), including outer vertical resonance (OVR) at u> = y/nO,± 
and inner vertical resonance (IVR) at Co = —^/nQ±. The VRs are turning points of eq. (29). The IVR exists only for m > sfn. 
Note that for Keplerian disks and n = 1, the LR and VR are degenerate. 
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For uj < 0, a Lindblad resonance (LR) exists only for m = 0, where u> = —k, and a vertical resonance (VR) may exist for 
m < y^n, but there is no corotation resonance in the disk. 

From the dispersion relation we can identify the wave propagation zones for ui > (see Fig. 1 and Fig. 2): (1) For n = 0, 
the wave zone lies outside the OLR and ILR (i.e., r > tolr and r < tilr); (2) For n > 2, the wave zones lie between ILR 
and OLR (tilr < r < tolr) and outside the IVR (if it exists) and OVR (r < tivr and r > rovn)\ (3) For n = 1, waves can 
propagate everywhere. 

The group velocity of the waves is given by 

C 9 = ~JT = —r T" ( 36 ) 

aK Co 1- ( K /Co) 2 (nQl/Cb 2 ) 

The relative sign of c g and and the phase velocity c p — co/k is important for our study of wave excitations in §§4-7: For u> > 0, 
positive (negative) c g /c p implies that trailing waves (k > 0) carry energy outward (inward), while leading waves (k < 0) carry 
energy inward (outward). The signs of c g /c p for different propagation regions arc shown in Figs. 1-2. Note that for n = and 
n > 2, c p — > oo and c g — > as k — > at the Linblad/ vertical resonances. But for n — 1, we have 

C9 = ^T^ sgn( ~ fc ^ ) - (37) 

Thus \c g \ — > c/2 at the Lindblad/vertical resonances §. 



4 LINDBLAD RESONANCES 

We now study wave excitations near a Lindblad resonance (where D = 0). Equation (29) shows that different n-modes are 
generally coupled. However, when solving the equation for a given n-mode, the coupling terms can be neglected if \rj n \ > \ rj n ±2 \ ■ 
Note that in the vicinity of a Lindblad resonance, \d\nD/dr\ S> \d\na/dr\ ~ \d\nQ./dr\ ~ 1/r. For a thin disk, if m and n not 
too large (m,n <§C r/h), the terms proportional to c~ 2 oc h~ 2 are the dominant non-singular terms. Keeping all the singular 
terms (oc D~ x ), we have 



d 2 /dlnD\ d 2mQ. fdlnD\ D(Cb 2 - nQ. 



d z /dlnD\ d + 2mf2 / 
dr 2 \ dr J dr rui V 



dr 

+ - — i — \nr) n + n + l)(n + 2 r/, l+2 = — — tp n ^~4>n, 38 

r \ dr J rdr c 2 lu 2 



where 



Now the terms proportional to (d\nD/dr)rj n and (dlnD/dr)7/ n+ 2 can be dropped relative to the other terms (see Goldreich 
& Tremaine 1978)^, we then obtain 



" d 2 


/dlnL>\ 


dr 2 


{ dr ) 



d D(u 2 - ntlj 
dr c 2 H) 2 



Vn = — j— tpu ^r<t>n- (40) 



We now proceed to solve eq. (40). Note that besides its (apparent) singular behavior, the resonance point of D — is a 
first-order turning (or transition) point when n^l, and a second-order turning point when n — 1 for Keplerian disks. These 
two cases should be investigated separately. 

4.1 n/1 Mode 

In the vicinity of the Lindblad resonance tl, we change the independent variable from r to 

x = (r — r L )/r L (41) 
and replace D by (r dD/dr) rL x + (r 2 d 2 D/dr 2 ) rL x 2 /2. For \x\ <S 1, eq. (40) becomes 

§ Of course, eqs. (35)- (36) are valid only away from the resonances, so the limiting values discussed here refer to the asymptotic limits 
as the resonances are approached. 

^ To see this explicitly, we set r? n = D 1 / 2 y„ and reduce eq. (38) to the form d 2 y„/dr 2 + f(r)y n = ■ ■ ■. We then find that the (d In D/dr)r] n 
term in eq. (38) gives rise to a term oc (din D /dr) oc D _1 in f(r), while the d 2 rj„/dr 2 and (din D/dr)dr) n /dr terms in eq. (38) both 
contribute terms proportional to (dlnD/dr) 2 oc D~ 2 in /(r). 
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Figure 1. A sketch of the function G = D(l — nD. 2 L /Cu 2 ) as a function of r for m = 2 and different values of n (all for u> > 0). The 
dispersion relation is G = —k 2 c 2 , and thus waves propagate in the regions with G < 0. The ± gives the sign of c g /c p of waves in the 
propagation region. 



1 d 



VJn 



- f3x)r) n = ax - jx , 

y dar x ax J x 

where ip n on the right-hand-side is evaluated at r — tl and 



(42) 



^dD/dr 



h 2 Cj 2 



{u 2 -nQ 2 ± )r 3 dD/dr 



7 = nr L 



j2 n / j„2 



h 2 n 2 ± L 



% d 2 D/dr 2 dD d 
2h 2 Cb 2 dr dr 



( K \ 



(43) 



Here we have kept the term —jx 2 in the Taylor expansion of the second term of right-hand side of Eq.(40) because the 
leading order term —ax generates only non-wave oscillations, as we shall see later. In orders of magnitudes, |/3| ~ (tl/K) 2 
and \a\ ~ (7! ~ n(rL/h) 2 \(f> n \. In particular, for disks with Q.± — k = Q we have 



/ 8 = 2(l-n)(l±m) 



r 2 dlnfi 
h 2 dlnr 



(44) 



where upper (lower) sign refers to the OLR (ILR). Note that ILR occurs only for m > 2, so the factor (1 ± rn) is never zero. 
To solve eq. (42), it is convenient to introduce a new variable fj defined by r\ n — drj/dx (see Ward 1986). We find 
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\ n=0 
CR \OI_R 


n=1 

CR OLR/OVR 


n=2 

CR OLRAOVR 


■ • • 

\ + 











Figure 2. Same as Fig. 1, except for m = 1. 



1 d 2 fi 
dx \ x dx 2 



13— = — - a--yx. 
dx x l 



Integrating once gives 



d 2 n . 
d^- pxr) 



-ipn — ax — -■yx + cx, 



where c is an integration constant. Then, let 
a 7 2 . c 

we have 

d 2 y „ 7 

As we see below, introducing the variable t/ singles out the wave part from rj n . The homogeneous form of Eq.(48) is the Airy 
equation and its two linearly independent solutions are Airy functions (Abramowiz & Stegun 1964, p. 446): 

yi = Ai(f3 1/3 x), y 2 = Bi(f3 1/3 x). (49) 

By the method of variation of parameters, the general solution to the inhomogeneous equation (48) can be written as 



(45) 

(46) 
(47) 

(48) 



// - ; 1/2 J .</: ^-'Ix - 2/1 



V2^dx + M Vl + Ny 2 , 



(50) 



where M and N are constants which should be determined by the boundary conditions, and W = yidy2/dx — y2dyi/dx 
/tt is the Wronskian. After writing £ = p 1 ^ 3 !, eq. (50) becomes 



V = 



£2/3 



r 



Ai&dt + N Bi(£) 



(51) 



where we have absorbed a factor /3 2 ^ 3 /(7rVl' n ) into the the constants M and N. 

The constants M, N in eq. (51) are determined by boundary conditions at |£| ^> 1. Note that the condition \x\ = 
|/?~ 1/3 £| < 1, which is the regime of validity for the solution described in this subsection, implies that |£| <C |/3| 1,/3 . Since 
\(3\ 3> 1, it is reasonable to consider the £| ^> 1 asymptotics of eq. (51). In the following we sometimes denote this asymptotics 
as |£| — > oo for convenience. 

For £ 3> 1, the asymptotic expressions of the Airy functions and their integrals are (Abramowiz & Stegun 1964) 



A 1 -l/2 e -l/4 -2e 3 / 2 

MO ~ 2 71 " ^ e 



Jo 



/IVi^Jt 1 1 -l/2 t -3/4 -|e 3 / 2 1 
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Lli(U = n-^ 2 r 1/4 e^ 3/2 - oo, f Bi{t)dt~ ^C^e^' 2 -> oo. (52) 

Jo 



Since grows exponentially when £ — > oo, the coefficient [J^ Ai(£)d£,+ N] before it in (51) must be very small, otherwise 

the quantity y (and hence r) n ) will become exponentially large, making it impossible to match with any reasonable physical 
boundary conditions. Without loss of generality, we will take this coefficient to be zero, based on the observation [see eq. (54) 
below] that the solution on the £ < side is nearly unaffected whether the coefficient is small or precisely zero. Thus 
N = — Ai(£)d£ = —1/3. Note that although J Q Bi(£)d£ in (51) also exponentially grows when £ — > oo, it is canceled by 
Ai(£) which is exponentially small. As the Airy functions are monotonic for £ > 0, wave solution does not exist on the side 
of £ > 0, or 

(1 — n)(l ± m)x < (Nonwave region) (53) 

[see eq. (44)]. This is consistent with the wave propagation diagram discussed in §3 (see Figs. 1-2). 

Now let us examine the £ < region. For £ <C — 1, the asymptotic behaviors of the Airy functions are 



I 



4 Ai(Od£ « -\ + jt-^-O- 3 / 4 cos A^), 

6 



Bi(£)^- 1/2 (-?r 1/4 cosAm 

f Bi(eK M - 1/2 (m 3/4 s mAm (54) 
Jo 
where 

A(0 = |(-C) 3/2 + J. (55) 
Equation (51) with N = -1/3 yields, for £ < -1, 

V - -^^(^r 1/4 [(! - i M y X(() + (1 + iM)e"««)] . (56) 
From the relation between r) n and y, we then obtain the asymptotic expression for r\ n at £ <g —1 as 

- f + ^ + iZ ^(^) 1/4 [(1 ~ ^)« <XK) - (1 + iM)e-«<«] . (57) 

The first two terms in eq. (57) describe the non-propagating oscillation, which is called non-wave part by GT. The last term 
gives traveling waves. Eq. (57) explicitly shows that waves exist on the side of £ = j3 1 ^ z x < 0, or (1 — n)(l ± m)x > 0. Again, 
this is consistent with the propagation diagram discussed in §3: If to > 0, the wave zones are located on the outer side of the 
OLR and inner side of the ILR for n = 0, whereas for n > 2, waves exist on the inner side of the OLR and outer side of the 
ILR. If lo < 0, LR is possible only for m = 0, and the wave zone lies in the outer side of the resonance for n = and appears 
on the inner side for n > 2. 

To determine the constant M, we require that waves excited by the external potential propagate away from the resonance. 
The direction of propagation of a real wave is specified by its group velocity, as given by eq. (36). For the waves going away 
from the resonance, we require sgn[c 9 ] = sgn[a;], and thus the local wave-number k must satisfy sgn[fc] = sgn[a;a;(l — n)]. On 
the other hand, the wavenumber associated with the wave term e ±tX( -^ in eq. (57) is k = ±dX/dr = =F/3 1,/3 (— £) 1 / 2 r~ 1 , which 
gives sgn[fc] = =psgn[/3] = ±sgn[a;] (since £ < 0). Accordingly, if ui > 0, for the n = OLR and the n > 2 ILR, the e lX term 
represents the outgoing wave, and in order for the e~ lX term to vanish, we have M = i, which leads to 



r?n^r? + 77=^ + ^a; + i % (-g) 1/4 e' X( °, (r > r OLR for n = or r > r ILR for n > 2) (58) 

where fj and fj represent the non-wave part and the wave part, respectively. Similarly, for the n = ILR and the n > 2 OLR, 
the e~ lX term represents the wave propagating away from the resonance, and eliminating the unwanted e lX term requires 
M = —i, which yields 

Vn^V + V=^ + ^3;- i 71 " (-C) 1/4 e~' x(0 , (r < r ILR for n = or r < r OLR for n > 2). (59) 

Lindblad resonance also occurs for u < and m — 0, in which case we find that eq. (58) applies for n > 2 and eq. (59) for 
n = 0. 

We can now use eq. (33) to calculate the angular momentum flux carried by the wave excited at the LR. Obviously, the 
m — mode carries no angular momentum. For the n = OLR and the n > 2 ILR, we substitute eq. (58) in eq. (33), and 
find 
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F„(r > r LR ) = — nlm-n' 



rdD I dr 



where, according to eqs. (39), (43) and (48), 



^ I d(f> n ^ 2mQ. 

I dr cu — mQ. 



+ n 



In 



2(k 2 -nQ 2 ± ) dr v /i 4 k 4 dr 



dD 



(60) 



(61) 



Similarly, for n = ILR and n > 2 OLR, we use eq. (59) in eq. (33) to find 



F„(r < r LR ) = n! run 



rdD I dr 



(62) 



Obviously, the angular momentum flux at the £ > side (where no wave can propagate) vanishes. Note that the torque on 
the disk through waves excited in r > r LR is F n (r > r LR ), while for waves excited in r < r LR the torque is — F„(r < r LR ). 
Since dD/dr < for OLR and > for ILR, we find that for both n = and n > 2, the total torque on the disk through both 
IRL and OLR is 



T„(OLR and ILR) = \F ORL | - \F ILR | = -n! rrni 2 



+ 



(63) 



That is, independent of n, the torque on the disk is always positive at OLR and negative at ILR (see GT). We note that for 
n = /j, = 0, our result agrees with that for the 2D non-self-gravitating disks (GT; Ward 1986). 



4.2 n — 1 Mode: Lindblad/Vertical Resonances 

For n = 1, the LR and VR are degenerate for a Keplerian disk, and we shall call them Linblad/ vertical resonances (L/VRs) 
II. The resonance radius r = r_L is both a (apparent) singular point and a second-order turning point of the wave equation 
(40). The wave propagation diagram (see Fig. 1-2) shows that waves exist on both sides of a L/VR. 
Expanding eq. (40) in the vicinity of the resonance (\x\ <g; 1), we have 

j2 



d Id ^,2 2 

dx* x dx 
where 



a\ = 



■ipi 



aix, 



r 3 ct>idD/dr 


, b = 


r 2 dD/dr 


h 2 £j 2 




hn 2 ± 



tpl = 



( d 2mil \ 

Kdr 5 Mj 0i - W; 



In orders of magnitudes, a?ij ~ (ri,/ft) 2 0i and b ~ Th/h. Substitution of 771 = y — tpix into eq. (64) gives 

d 2 Id 2 2 ,2,3 

— — y — y + b x y = b tpix - aix. 

dx 2 x dx 

The two linearly independent solutions to the corresponding homogeneous form of eq. (66) are 

-ibx 2 /2 ibx 2 /2 

Vi = e , 2/2 = e 

The method of variation of parameters then gives the general solution of eq. (66) : 



y = e 



:C 2 /2 



iC 2 /2 



S(0dC + N 



+ e 



i ( 2 /2 



< 2 /2 



5(C)dC + M 



(64) 



(65) 



(66) 



(67) 



(68) 



where we have defined £ = h 1 * 12 ! and 



Note that although our analysis is limited to \x\ <S 1, we have extended the integration limit to £ = ±00 in eq. (68) — This 
is valid because t~ri//!> 1 for a thin disk and the integrands in the integrals are highly oscillatory for |£| ^> 1 (so that the 



Bate et al. (2002) previously analysed the mixed L/VRs for axisymmetric waves (m = 0); such waves do not carry angular momentum. 
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contribution to the integrals from the |£| 3> 1 region is negligible; see Wong 2001). The wave solution in the |£| 3> 1 region is 
approximately given by 



V 



y = 



-i C 2 /2 



S(t)dC + N 



+ Me 



-iC/2 



« » 1) 



S(t)dC + M 



-«'/*+ Ne**' 2 , (C«-l). 



(70) 



(71) 



The constants M and N can be fixed by requiring that no waves propagate into the resonance. From our analysis of the wave 
group velocity in §3 (see Figs. 1-2), we find that for the waves propagating away from the resonance, the local wavenumber 
must be positive, irrespective of whether it is inner or outer n = 1 L/VR. Accordingly, we must have M — [from eq. (70)] 
and N = [from eq. (71)]. Since the integral 



/ 



ir/4 



cqs. (70)-(71) reduce to 

v " \/2i-^ +i m e2 4 ' (c>>1) ' 



y ~ V 2 V6V2 + *&3/2 

Using eqs. (33) and (73), we find that the angular momentum flux carried outward from the resonance (toward larger 
radii) by the wave excited at a L/VR is given by 



-<3C 2 +»f 



(C « -i) 



(72) 

(73) 
(74) 



ra<f>l 
h 2 dD/dr 



(75) 



where in the second equality we have used the fact |ai/b| = (r/h)\<f>i\ ^> |^i|. Similarly, using (74), we find that the angular 
momentum flux carried by the x < wave is 



p i( r < r L,vn) = 



(. 



ra<j>'( 



L/VR > 2 \h 2 dD/dr 



(76) 



Thus the angular momentum transfer to the disk through a L/VR is simply Ti(L/VR) = Fi(r > r L/VR ) — F\{r < r L/VR ) = 
2Fi(r > r L/VR ). Combining the inner and outer resonances, the total torque is 



Ti(OL/VR and IL/VR) 



( //'-/ 



/d V OL /VR WdD/dr) 



IL/VR 



(77) 



Again, we find that the torque on the disk is positive at OL/VR and negative at IL/VR. Comparing the above result with 
the results of §4.1 and §5, we find that although the n = 1 L/VR involves a combination of LR and VR its behavior is more 
like a VR. 



5 VERTICAL RESONANCES (TV > 2) 

We have already studied the n = 1 vertical resonance in §4.2. We now examine VRs, Co 2 = rani, for ra > 2. 

In the neighborhood of a VR, there is no singularity in eq. (29). For a thin disk (h <JC r) with m,n <§C (r/h), it is only 
necessary to keep d 2 rj„/dr 2 and the terms that are oc c~ 2 oc h- 2 . This can be justified rigorously from the asymptotic theory 
of differential equations (e.g., Olver 1974), which shows that the discarded terms have no contribution to the leading order 
of the asymptotic solution. Indeed, the discarded terms only make a small shift to the vertical resonance for a thin disk. As 
for the coupling terms with other modes, in addition to the reasons given in §4, the dropping of the coupling terms are more 
strongly justified here as the VRs with different ra's are located at different radii and their mutual effects can be considered 
insignificant. Therefore, around the VR radius r v , we can simplify eq. (29) to 

d 2 D{Cb 2 -ntij) nD 

dr^ n h^n\ h^ n - ( ™> 

Changing the variable from r to x = (r — r v )/r v , we obtain, for |a;| <§C 1, 
d 2 

j^Vn - Aa;?7n = X, (79) 
where 
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2r 3 D(mCjdQ/dr + nSl±dSl±/dr) 



x = — 



nr D(j> n 



h 2 u 



Similar to §4.1, the general solution to eq. (79) reads 



^X 

A 2 / 3 



Ai^)(k + N) Bi{<;) 



-a 



Bi{<;)(k + M) Ai{s) 



(80) 



(81) 



where ? = \zx. Suppressing the mode which is exponentially large when ? — » co yields N = —1/3. 

Waves can propagate in the ? < region, i.e., on the outer side of the OVR and on the inner side of the IVR (see 
Figs. 1-2). For the waves to propagate away from the resonance, we require the group velocity to satisfy sgn[c 9 ] = sgn[x], and 
hence the local wavenumber to satisfy sgn[fc] = sgn[a;ai(l — n -1 )]. Thus, if u > 0, we demand that as <; 



for the OVR and r) n oc e 



i(-0 3/2 



GO , Tj n OC 6 3 q > 



for the IVR. This determines M and yields, for ? 



exp 



-i(^(-.) 3/2 + T) 



The angular momentum flux is then 

,2 



F n (r > r r 



n\ run 



F„(r < r IVR ) = -n\mn 



°X 
\D 

vx~ 
\D 



n\\fn 



( r > r ovn.) 
( r < r IVR ). 

(- 



y/7i + m \ti 2 £ldQ,/dr 

2 



OVR 



raq 



— 72 ! V 72 

2 y/ri — m \h 2 QdQ,/dr 



(82) 
(83) 

(84) 
(85) 



The torque on the disk is F n (r > r OVR ) at OVR and — F n (r < r IVR ) at IVR. Note that the IVR exists only for m > \Jn (see 
§3), so — F n (r < r IVR ) < 0. The total torque on the disk due to both OVR and IVR is 



T n (OVR and IVR) = -— 



s/n + rn \h 2 QdQ/dr 



ra<p„ 



(86) 



ovr \fn — m \h 2 SldQ,/dr ^ 

(Obviously, in the above expression, the IVR contribution should be set to zero if m < y/n). Again, we see that the torque is 
positive at OVR and negative at IVR. 

If w < 0, a single VR exists for m < \fn. In this case waves are generated on the outer side of the resonance with 



Vn(r >r v )- 
as c = A 1/3 x 



* 1/2 x 



A2/3 

-co. Thus 



(-?) * exp 



-(|(-,f /2 + J) 



(r > r v \ for u> < 0) 



F n (r > r v ) = —n\y/n 



(- 



yfn — m \h 2 Q,d£l/dr 



(w < 0). 



(87) 



(88) 



The torque on the disk due to such VR is negative. 

It is interesting to compare our result with the one derived from shearing sheet model (Takeuchi & Miyama 1998; Tanaka 
et al 2002): 



F n = ^nU 



h 2 n\dn/dr\ 



(89) 



Clearly, our eq. (84) reduces to eq. (89) for m >• \Jn. 

At this point, it is useful to compare the amplitudes of the waves generated at different resonances. For LRs, F n ~ cr^/fi 2 ; 
for both the n = 1 L/VRs and n > 2 VRs, F n ~ (r /h) 2 acj>^/Q 2 . Thus, when the external potential components have the 
same orders of magnitude, the angular momentum transfer through VRs is larger than LRs for thin disks. Since we expect 
4> n oc (h/r) n , the n = 1 vertical resonance may be comparable to the n = Lindblad resonace in transferring angular 
momentum. 



6 COROTATION RESONANCES 

Corotation resonances (CRs), where u = or to — mfi, may exist in disks for lo > 0. The WKB dispersion relation and 
wave propagation diagram discussed in §3 (see Fig. 1-2) show that for n = 0, waves are evanescent in the region around the 
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corotation radius r c , while for n > wave propagation is possible around r c . This qualitative difference is also reflected in the 
behavior of the singularity in eq. (29). We treat the n = and n > cases separately. 



6.1 n = 0Mode 



In the vicinity of corotation, we only need to keep the terms in eq. (29) that contain the uj~ 1 singularity and the term oc h~ 2 . 
The term oc r)n+2/£j can also be dropped, since from eq. (26) we can see that the coupling term is negligible when \Q\ is small. 
Thus for n — 0, eq. (29) can be approximately simplified to 

j2 



d 2 2mft / d 

-no 



dr 2 ' 1 " rCu Wr D J m h 2 tl 2 ± ' 
Near the corotation radius r c , we introduce the variable x = (r — r c )/r c , and eq. (90) becomes 
d 2 r] 



D 



-no 



2mQ ( d , aQ,\ , 4umfi , 

- — ( ~r m tt ) ^° + -^r-fo- 

tuj \dr D J r z uj 



dx 2 
where 



no 



x + it 



P = 
$ = 



20 d , /o-n\ 



dQ./dr dr 
<t>o + 



In 



\ D J 



Dr 2 



h 2 Q, 2 



1/2 



= (Wc) rc > 1, 



rdln(crn/D)/dr 



(90) 



(91) 



(92) 



(93) 



In eq. (91), the small imaginary part it (with e > 0) in \/x arises because we consider the response of the disk to a slowly 
increasing perturbation (GT) or because a real disk always has some dissipation which has not been explicitly included in our 
analysis. If fi = or |0 2 | <C \4>o\, then the eq. (91) is the same as the equation studied by GT for 2D disks. 

Goldreich & Tremaine (1979) solved eq. (91) neglecting the (p/x)r)o term, and found that there is a net flux of angular 
momentum into the resonance, given by 



AF C = F(r c -) - F(r c +) = 2tt m 



<1> 2 



dQ/dr dr\~D) 



(94) 



On the other hand, integrating eq. (91) from x = 0— to x — 0+ gives (drio/dx)o+ — (drio / dx)o~ = iirp(rio + $), which, together 
with eq. (32), yield (Tanaka et al. 2002) 



AF C 



Wo 



$M d foQ, 



dQ, I dr dr 



crS2\ 

~D ) 



(95) 



Tanaka et al. (2002) argued that neglecting rjo in eq. (95) may be not justified in a gaseous disk and that the revised formula 
(95) fits their numerical result better than eq. (94). 

Although eq. (95) is a precise result of eq. (91), the presence of the unknown quantity rjo(r c ) makes the expression not 
directly useful. Therefore, a more rigorous solution of eq. (91) without dropping the (p/x)rjo term seems desirable. In the 
following we provide such a solution. 

After changing variable w = rjo + <E>, eq. (91) becomes 

(96) 



TT + ( q )w = -q <S> 



dx 2 

To solve this non- homogeneous equation, we need firstly to find the solutions of the corresponding homogeneous equation 

By introducing the variables 
w — xe qx y, s — —2qx, 

eq. (97) is transformed to Rummer's equation (Abramowitz & Stegun 1964) 

0. 



(98) 



Its two independent solutions can be chosen as 



yi = U(l + ^,2,s), V 2 = e s U{l 



2q 



,2,-s), 



(99) 



(100) 



© 2005 RAS, MNRAS 000, 1-20 



14 Hang Zhang and Dong Lai 



where the function U is a logarithmic solution denned as 

(_!)«+! 



U(a, n + 1, z) = 



nir(o-n) 



M(o, n + l,z) In 2 



E 



(a) r 2 r 
(n + l) r r! 



{i/>(a + r) - V(l + r) - + n + r)} 



— 

(n ~ 1)' -n ( a _ rl ) 
■ T(a) 



\ ^ (a 
^ (1 



(l-n) r r\ ' 



(101) 



in which M(a,b,z) is Rummer's function, i/> is the Digamma function and (...)r denotes the Pochhammer symbol (see 
Abramowitz & Stegun 1964). So the two independent solutions to eq. (97) are 



(102) 



(103) 



«Ji = xe qx U(l + |-, 2, -2ga;), w 2 = xc-«"l7(l - ^, 2, 2^), 
and their Wronskian can be shown to be 

d d 1 i,r(l-P) 

da; da; 2g 

Note that owing to the singularity of the function U at the branch point x = 0, a branch cut has been chosen in the lower 
half of the complex plane. 

The solution to eq. (96) has the form 



W2 



/x 
-o 



wi- 



da; + wi 



W 

= 2q i ^e- t7 ' {1 -^ ) xe- qx U(l~ £-,2,2qx) [ xe qx U(l + 2, ~2qx)dx 



W2 — TT7—dX 



w 



+2q*<S>e- 17r(1 -% ) xe qx U(l + £-,2,-2qx) / x e - qx U(l - 2, 2qx)dx, 

2q J x 2q 



(104) 



where we have adopted the physical boundary conditions that for \qx\ 3> 1, the quantity |w| does not become exponentially 
large. One of the limits of each integration has been extend to infinity since it only introduces negligible error due to the 
property of the integrand. From eq. (101), at x = 0, we have 



w(x = 0) = e 

4 



e~^tU(l + £-,2,i)dt 
2q 



2q> Jo 

4 6 r(l + ^)X W{1 2q^^ )dt 



(105) 



Utilizing a formula of Laplace transform (Erdelyi et al. 1953) 



f 

Jo 



e - st t"- L U(a,c,t)dt : 



r(6)r(6- c + i) 



T(a + b-c+l) 
to calculate out the integrals, we obtain 

pn 



2 F 1 (b, b - c + 1; a + b - c + 1; 1 - s) 



w(x = 0) = — — e 1 



2, ' 



2q 



7T 



Fi(2,l;2+£;i) 



2^1(2,1; 2- 



2q ' 2 / 



1 + £ - 

^ 2q 



1 



2q 



where 2.F1 denotes the Gaussian hypergeometric function. Hence, 



1 £I Sm 2? 

Vo(r c ) = -e -g— 



and from eq. (95), 



i*i(2,l;2 + £;i) 3^(2,1:2- 



1 + -E- 



1 _ _E_ 

2q 



(106) 



(107) 



(108) 



AF C = liCm 



$ 2 d (oQ, 



dQ/dr dr l~D~/ 



16 (p. 



I S in 2 | r 2 F 1 (2,l;2+^;i) +!J F 1 (2,l;2-^;i)^ 



1 + 



1 - -E- 

2 9 



(109) 



By the fact 2-Fi(2, 1; 2; i) = 2, it is easy to check that when p/q ~ /i/r — > 0, |/?o(tc)| becomes much smaller than |$|, and AF C 
reduces to eq. (94), the original Goldreich-Tremaine result. 
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6.2 n > 1 Mode 

In the vicinity of corotation with n > 1, the terms oc h~ 2 in eq. (29) are dominant, and we only need to keep these terms and 
the second-order differential term. Eq. (29) reduces to 

d 2 D(Q 2 -nQl) nD 

where <f> n is evaluated at r = r c . Defining x = (r — r c )/r c and expanding eq. (1 10) around x = 0, we have 

d ^ + C-^— i r ]n = ^C T ^— 1 (111) 
dx 1 (x + ie) 2 (x + ie) 2 

where 

n 

2 



c= A 



hdQ/dr 



We remark here that eq. (110) is similar to the equation derived in the context of the stability analysis of stratified flows 
(Booker & Bretherton 1967), although the physics here is quite different. 
The general solution to eq. (Ill) is 

Vn = -4>u + Mz 1/2 z lv + Nz 1/2 z-™ = -<f> n + Mz 1/2 e l " lnz + Nz 1/2 e~ wlnz , (113) 

where v = UC — \ 1, z = x + ie (with e > 0) and M and N are constants. The first term, the non-wave part, is a particular 
solution, while the other two terms are solutions to the homogeneous equation, depicting the waves. 

From eq. (36), we find that the group velocity of a wave (with wave number k) near corotation is given by c g — 
—CLb 2 /(^/nKQ,±_)sgn(kCb). Thus sgn(c 9 ) = — sgn(fca;) for \x\ > 0. The z x l 2 z %v component in eq. (113) has local wave number 
k = d(vlnz)/dr = vj(r c x), and thus it has group velocity c g < 0. If we require no wave to propagate into x = from the 
x > region (as we did in studying LRs and VRs; see §4 and §5), then we must have M = 0. Similarly, requiring no wave 
to propagate into corotation from the x < region gives iV = 0. Thus we have shown explicitly that waves are not excited 
at corotation. Further calculation shows that even adding the higher-order terms to eq. (Ill) does not alter this conclusion. 
This is understandable because |fe| — > oo near corotation, and short-wavelength perturbations couple weakly to the external 
potential. 

However, unlike the n = case, waves with n > 1 can propagate into the corotation region and get absorbed at the 
corotation. To calculate this absorption explicitly, let us consider an incident wave propagating from the x > region toward 
x = 0: 



c 1 / 2 e ll/lnx , (x > 0; incident wave), (114) 



where A+ is constant specifying the wave amplitude. To determine the transmitted wave, we note that z = is the branch 
point of the function e lvlnz , and the physics demands that the solution to eq. (Ill) be analytic in the complex plane above the 
real axis. Thus we must decide the branch cut of the function so that it is single- valued. As discussed before (see §6.1), while 
our analysis in this paper does not explicitly include dissipation, a real disk certainly will have some viscosity, and viscous 
effect can be mimicked in our analysis by adding a small imaginary part ie (with e > 0) to the frequency lo. Alternatively, 
this can be understood from causality where the disturbing potential is assumed to be turned on slowly in the distant past. 
This is the origin of imaginary part of z — x + ie in eq. (Ill) or eq. (113). Therefore, we can choose the negative imaginary 
axis as the branch cut of the function e ll,lnz . In doing so, we have e ll,lnz = e lI/( - ln l a: l+ , ' r ) for x < 0. Thus the transmitted wave 
is given by 

rj n = i A+ e^ 7rv (— x) 1 ^ 2 e* I/ln ^ _a! \ (x < 0; transmitted wave). (115) 

Since v 3> 1, the wave amplitude is vastly decreased by a factor e~™ after propagating through the corotation. From eq. (32) 
or (33), the net angular momentum flux absorbed at corotation is 



AF C = F n (r c -) - F n (r c +) = n\nm (^) v\A+\ 2 (l + e' 2 ™) 



^^H^^pwJ^ 1 ' (116) 

where in the last equality we have used v ~ \/~C ^> 1. 

Similarly, consider a wave propagating from x < toward x — 0: 

r h = A-z 1/2 e- lvlnz =iA-e lxv {-x) 1/2 e- ivl ^- x \ (x < 0; incident wave). (117) 

The transmitted wave is simply 
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rj n = A-X 1/2 e~ il,lnx , (x > 0; transmitted wave). (118) 
The net angular momentum flux into the corotation is 

AF C = -„! ™ {^jlA-e-fd + e — ) * -n! ^ (^^) ^ e"f, (119) 

where the negative value of AF C arises because waves inside corotation carry negative angular momentum. 

In summary, for n > 1, waves propagating across the corotation are attenuated (in amplitude) by a factor e~ nl/ . Thus 
the corotation can be considered as a sink for waves with n > 1 — a similar conclusion was reached by Kato (2003) and Li et 
al. (2003), who were concerned with the stability of oscillation modes in accretion disks around black hole. The parameters 
A+, A- in eqs. (114) and (117) are determined by boundary conditions. In Appendix B, we discuss a specific example where 
waves excited at Lindblad/ vertical resonances propagate into the corotation and get absorbed and transfer their angular 
momenta there. 



7 WAVE EXCITATIONS AT DISK BOUNDARIES 



In the preceding sections we have examined the effects of various resonances in a disk. Even without any resonance, den- 
sity/bending waves may be excited at disk boundaries. 

To give a specific example, let us consider the n = m = 1 bending wave in a Keplerian disk driven by a potential with 
frequency co in the range < w < fl(r out ). Here we use r out and r in to denote the outer and inner radii of the disk. This 
situation occurs, for example, when we consider perturbations of the circumstellar disk of the primary star driven by the 
secondary in a binary system. Since no resonance condition is satisfied for the whole disk, the general solution to the disk 
perturbation is given by [see eq. (A9) in Appendix A]** 



r?i =Q- 1 G+(D/ra) 1/2 Q 
where 



-1/4 



Mexp(— i 



f Q 1/2 



dr) + N exp(i 



/ Q 1/2 



dr) 



Q = k 2 = 



h 2 Cb 2 Q. 2 



G = - 



D 



(120) 



(121) 



To determine the constant M, N, we assume that the inner boundary is non-reflective, thus N = 0. For the outer boundary, 
we assume that the pressure perturbation vanishes, i.e., 771 = 0, this determines M and eq. (A9) then becomes 



m = q~ 1 g- 



Q-^G 



(D/ra) 1/2 Q- 1/4 exp(i 



Q 1/2 dr). 



(D/ray/zQ- 1 / 4 

A direct calculation using eq. (33) shows that the angular momentum flux carried by the wave is 
roG 2 



Fi = TV 



Q' i / 2 D 



raD<f> 2 
hWQ 3 / 2 



(122) 



(123) 



rout x 7 runt 

Therefore, in this model, waves are mainly generated at the outer disk boundary, propagating inward, while the angular 
momentum is transferred outward. 



8 DISCUSSION 

Here we discuss several assumptions/issues related to our theory and how they might affect the results of the paper. 

(i) Radially Nonisothermal Disks. In deriving the basic equation (29) for the disk perturbation, we have dropped several 
terms proportional to the radial gradient of the sound speed (see footnotes 1 and 2). It is easy to see that these terms vary 
on the lengthscale of r and do not introduce any singularity or turning point in our equation, therefore they do not affect 
any of our results concerning wave excitation/damping studied in §§4-7. However, a r-dependent sound speed gives rise to a 
nonzero dfl/dz, which can modify the structure of the perturbation equation near corotation. Indeed, if the (unperturbed) 
surface density a and sound speed c profiles satisfy simple power-laws, a oc r~ a and c oc r^ 13 , then the angular velocity profile 
in the disk is given by (Tanaka et al. 2002) 



** The wavenumber k for the n = m = 1 mode (for Keplerian disks) is given by k 2 h 2 = (u/n) 2 (2n - w) 2 /(n - w) 2 , which reduces to 
k 2 h 2 = (2uj/U) 2 -C 1 for u) < Q(r ou t) <C CI. Thus the radial wavelength may be much larger than h and comparable to r, in which case 
the WKB solution is not valid. 
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fi = Q.K < 1 



(124) 



where Qk (r) is the Keplerian angular velocity. For a thin disk, the deviation of fi(r, z) from fix is obviously very small. 
Nevertheless, if the 2-dependence of fi is taken into account, an additional term, (—f3mnQ,D/r 2 u> 3 )w n should be be added to 
the left-hand-side of eq. (29)11. Obviously, this term does not affect waves near a Lindblad resonance and vertical resonance. 
Because of the strong u)~ 3 singularity at corotation (where w = 0), one might suspect that our result on the n > 1 corotation 
resonance (see §6.2) may be affected (see Tanaka et al. 2002). In fact, we can show this is not the case. With the iD~ 3 term 
included, the perturbation equation near a n > 1 corotation resonance is modified from eq. (Ill) to 

d 2 r) n C Ci C4> n Ci4> n 

+ ~>Vn + -tVn = -= (125) 



dx 2 x 2 x 3 x 2 x 

where 

n 2 n 



r 3 (dQ/dr) 3 



(126) 



Since C ~ (n/m 2 )(r c /h) 2 3> 1 while C\ ~ nf3/m 2 , the new terms are important only for \x\ £ /3(h/r c ) 2 . Thus we expect that 
our solution, eq. (113), remains valid for |x| 3> /3(h/r c ) 2 . Indeed, the general solution of eq. (125) is 

Vn = -K + Mx?J l2l/ (2C?x~ 5 ) + Nx? J- av (2C\ aT * ), (127) 



where v = yj C — 1/4, and J±2i^ is the Bessel function (Abramowitz & Stegun 1964, p. 358). This solution approaches the 
form of eq. (113) for \x\ ^> f3(h/r c ) 2 . Therefore, the analysis in §6.2 remains valid and our result on wave absorption at the 
n > 1 corotation resonance is unchanged. We conclude that while our theory is explicitly based on radially isothermal disks, 
our results remain valid when this condition breaks down. 

(ii) Vertical structure of disks. Our theory is concerned with vertically isothermal disks, for which 3-dimensional pertur- 
bations can be decomposed into various Hermite components [eq. (15); see Tanaka et al. 2002]. It would be interesting to 
investigate if similar decomposition (with different basis functions) is possible for more general disk temperature profiles. To 
simplify our equations, we have also neglected stratification in our analysis [see eq. (8)]. In particular, vertical stratification 
gives rise to a local Brunt- Vaisala frequency of order fix (z/h) 1 ^ 2 , and it is not clear to what extent such stratification will 
affect our results involving vertical fluid motion in the disk (see Lubow & Ogilvie 1998). This issue requires further study. 

(iii) N on- Keplerian Disks. Although in this paper we have considered Keplerian disks (for which fi = k — fix ) , extension 
to more general disks is not difficult. Indeed, since we have been careful to distinguish fi, k, fix throughout the paper, most 
of our equations are valid when fi ^ k ^ fix- The only exception is the n = 1 Lindblad/vertical resonances studied in §4.2: 
Only for a Keplerian disk (fix = k) is the n = 1 vertical resonance degenerate with the Lindblad resonance, and such a 
combined Lindblad/vertical resonance needs special treatment. For a disk where the Lindblad resonance and n = 1 vertical 
resonance are well separated, they must be treated separately, with the procedure similar to those given in §4.1 (for Lindblad 
resonances) or §5 (for vertical resonances). For a nearly Keplerian disk, with |fix — k\/k <§C 1, the LR (at tl) and the n = 1 
VR (at rv) are rather close, and the problem requires some consideration. Expanding eq. (40) around the Lindblad resonance, 
we find [cf. eq. (64)] 

d2 1 d , 1,2 i \ V>i ,- 10Q N 

-r-zVi j-Vi + o xix - xvmi = otix, (128) 

dx 2 x dx x 

where xy = (rv — rL)/rL and on, b, ipi are given by eq. (65). Obviously, solution (67) breaks down for \x\ <; \xv\- For 
bx v <^ 1, or \xv\ <C b~ x / 2 ~ (h/ri,) 1 ^ 2 , the asymptotic solutions (73) and (74) remain valid and the angular momentum flux 
is unchanged. For bx v <; 1, the angular momentum flux expression has the same form as eq. (75) or (76) except that the 
pre- factor (miv 2 /2) may be changed by a factor of order unity. 

(iv) Nonlinear effect. From the wave solutions we have derived in this paper, we see that the enthalpy or density per- 
turbation of the disk is finite at various resonances. However, the azimuthal velocity perturbation can become singular at 
the corotation resonance (GT). Viscous and nonlinear effects may become important at the resonance and affect the derived 
torque formula. Therefore, linear, inviscid theory for the corotation resonance is incomplete even for a very small external 
potential. As has been pointed out by Balmforth & Korycansky (2001), and largely discussed in hydrodynamic stability theory 
(see Drazin & Reid 1981), critical layer may emerge at the resonance. This issue requires further study (see Ogilvie & Lubow 
2003; Goldreich & Sari 2003). 



tt Nonzero dfl/dz also gives rise to an additional "coupling" term, oc (f3/oj 3 )ri n +2- This can be neglected [see the discussion following 
eq. (29)]. 
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9 CONCLUSION 

In the paper we have studied the linear response of a 3D gaseous disk to a rigidly rotating external potential. The disk 
is assumed to be non-self-gravitating, geometrically thin and vertically isothermal. The external potential and the disk 
perturbation can be decomposed into various Fourier- Hermite components, each proportional to H n (z/h) exp(im6> — iuit)], 
characterized by the azimuthal index m and the vertical index n which specifies the number of nodes along the 2-direction 
in the fluid density perturbation. We have derived analytical expressions for the various wave modes excited at Lindblad 
resonances and vertical resonances, and calculated the angular momentum fluxes associated with these waves and hence the 
torque on the disk from the external potential. We have also studied wave damping and angular momentum transfer at 
corotation resonances. For wave modes that involves only 2D motion (n = 0), our general formulae reduce to the standard 
result of Goldreich & Tremaine (1979). 

Our main results on wave excitation/damping can be most conveniently summarized using the wave propagation diagram 
(Figs. 1-2) which follows from the dispersion relation [eq. (35)]. In 2D disks, waves are excited only at the inner and outer 
Lindblad resonances, and propagate away from the corotation. By contrast, in 3D disks, additional channels of wave generation 
open up through vertical resonances, and waves can propagate into corotation and where angular momenta are deposited. 
Irrespective of the direction of propagation of the excited waves, the torque on the disk is positive for waves generated at 
outer Linblad or vertical resonances and negative at inner resonances. 

Our paper contains a number of analytical results which are obtained for the first time. A guide to the key equations 
are as follows: (i) Lindblad resonances: For the n — and n > 2 modes, the wave amplitudes excited at the resonances are 
given by eqs. (58)-(59), the associated angular momentum fluxes eqs. (60) and (62) and the torque on the disk eq. (63); for 
n=l, the Lindblad resonance and vertical resonance coincide for a Keplerian disk, the corresponding equations are (73)-(74), 
(75)-(76) and (77). (ii) Vertical resonances: For n > 2, the waves amplitudes excited at the vertical resonances are given by 
eqs. (82)-(83), the angular momentum fluxes eqs. (84)-(85) and the torque on the disk eq. (86). (iii) Corotation resonances: 
For n = 0, waves cannot propagate around the corotation region, but a torque is delivered to the disk at corotation. An 
improved expression for the torque is given by eq. (109), which reduces to the standard Goldreich- Tremaine result (94) in 
the h/r — > limit. For n > 0, waves can propagate into the corotation. The angular momentum flux deposited to the disk at 
corotation is given by eq. (116) or (119), depending on the incident wave amplitude. 

The last paragraph refers to waves excited by a prograde-rotating potential (with pattern speed fi p = ui/m > 0). It is of 
interest to note that for m < \Jn, vertical resonant excitations exist for a retrograde-rotating potential (with Q p < 0, i.e., the 
perturber rotates in the direction opposite to the disk rotation): The excited wave has an amplitude given by eq. (87) and 
angular momentum flux given by eq. (88). 

Even without any resonance, waves can be excited at disk boundaries. An example is discussed in §7. 

An interesting finding of our paper is that for a given potential component <j) n with n > 1, vertical resonances are much 
more efficient [by a factor of order (r/h) 2 ] than Lindblad resonances in transferring angular momentum to the disk. Whether 
vertical resonances can compete with n = Lindblad resonances depend on the relative values of (j> n (n > 1) and 0o- Since 
we expect (j>„ oc (h/r) n , the angular momentum transfer through the n — 1 vertical resonance may be as important as the 
n — Lindblad resonace when considering the perturbation of a circumstellar disk by a planet in an inclined orbit. We plan 
to investigate this and other related issues discussed in §8 in the future. 

APPENDIX A: SOLUTIONS AWAY FROM RESONANCES 

In the region away from any of the resonances, we only need to retain in eq. (29) the first and second derivative terms as well 
as the terms oc hT 2 . Thus for n > 0, eq. (29) reduces to 




(Al) 



where for n > 0, 




(A2) 



(A3) 



Q 



D 



(A4) 



h 2 n\' 
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2u / d d , D 2mfi\ , 
+ —[- r -- r ln- + —^)ct> 2 , n = 0. 

r \dr dr a no J 



(A5) 



Clearly, Q = k 2 gives the WKB dispersion relation (35), and \Q\ ~ h 2 . After making a transformation of the variable, 

y = (m/D) 1 ^, (A6) 

and using \Q\ 3> 1/r 2 , eq. (Al) reduces to 

pi + Qy = (ra/D)^G. (A7) 

In the wave-propagating region (Q > 0), the Liouville- Green or WKB approximation gives two independent solutions to 
the homogeneous equation of (A7) (e.g., Olver 1974) 



/r rr 
Q 1/2 dr), y 2 = Q~ 1/4 exp(i / Q 1/2 



dr). 



(A8) 



A particular solution to eq. (A7) can be written asymptotically, owing to the large parameter Q, as y = Q 1 (ra/D) 1 ^ 2 G, 
Thus the general solution to eq. (Al) is 



Vn = Q- 1 G+(D/ra) 1/2 Q~ 1/4 



/r pr 
Q 1/2 dr) + iVexp(i / Q 1/2 dr) 



(A9) 



where the two constants M and N are determined by boundary conditions. The boundary conditions may be imposed at the 
true disk boundary, or by the requirement to match the wave coming from a resonance. 

The the evanescent region (Q < 0), the general solution to eq. (Al) can be obtained in a similar way: 



r, n =Q- 1 G+{D/raf' 2 \Q\- 



1/4 



/r pr 
\Q\ 1/2 dr)+Nexp{ / \Q\ 1/2 dr) 



(A10) 



APPENDIX B: AN EXAMPLE OF WAVE EXCITATION AND ABSORPTION 



Here we give an example of the global solution of wave excitation at a Lindblad/vertical resonance, wave propagation, and 
wave absorption at corotation. For definiteness, we consider the n = 1 case, and assume that the wave excited at the L/VR 
(see §4.2) does not suffer any damping as it propagates toward the corotation radius. The wave excited from the outer L/VR 
at r — r OL is given by eq. (74) and can be rewritten as 



. /7T iir/4 

Vi = -ly/ 2 e 



h 1 / 2 \dD/dr\^ 



exp 



.1 


dD/dr 


~ l 2 


hQ 2 ± 



(r-r OL Y 



(Bl) 



For the region between the corotation and the outer L/VR, r c < r < r OL , the general solution is given by eq. (A9). Note that in 
this region, the group velocity c g has opposite sign as the phase velocity (see Figs. 1-2). Keeping only the inward- propagating 
wave term, we have 



rn =N{D/ray ,2 Q- L/ *exp{i I Q L/2 dr), 

where ro is a fiducial radius (r c < ro < r OL ). Matching (B2) with (Bl) as r — > r OL gives 



N 



I" ivr/4 

2' 



(, 



1/2 



h 1 6 1 



exp(- 



lV2, 



\\dD/dr\ 

Matching (B2) with r?i = A + x 1/2 e ivlnx [see eq. (114)] as r — > r c [or x = (r - r c )/r c -> 0+] yields 



/ hn\dn/dr\ 



2 ^ 



1/2 



(ro 
\dD/dr\ 



1/2 



h 1 (pi 



The angular momentum flux associated the wave at r = r c + is then [see eq. (116)] 



F 1 (r = r c +) 



Mi\dn/dr\ 



( ra<t> 



2 \h 2 \dD/dr\J r 



(B2) 



(B3) 



(B4) 



(B5) 



That is, the angular momentum absorbed at r = r c + is —Fi(r = r c +). Not surprisingly, this is exactly equal to the inward 
angular momentum flux transferred away from the OL/VR, i.e., —F\(r < r OL/VR ) [see eq. (76)]. 

We can similarly show that the angular momentum transferred by the waves excited at the inner Lindblad resonance and 
propagating to the corotation is absorbed at corotation, as well as the waves themselves. 
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